Prom incommensurate to dispersive spin-fluct nations: The high-energy inelastic 

spectrum in superconducting YBa2Cu306.5 
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We have investigated the spin fluctuations at energy transfers up to ~ 110 meV, well above 
the resonance energy (33 meV) in the YBa2Cu306.5 ortho-II superconductor using neutron time- 
of-flight and triple-axis techniques. The spectrum at high energies differs from the low-energy 
incommensurate modulations previously reported where the incommensurate wave vector is largely 
independent of energy. Well above the resonance the peak of the spin response lies at wave vectors 
that increase with energy. Within error the excitations at all energies above the resonance are 
best described by a ring around the (tt, tt) position. The isotropic wave-vector pattern differs from 
a recently reported square pattern in different but related systems. The spin spectral weight at 
high-energies is similar to that in the insulator but the characteristic velocity is ~ 40% lower. We 
introduce a method of extracting the acoustic and optic weights at all energies from time-of-flight 
data. We find that the optic spectral weight extends to surprisingly low-energies of ~ 25 meV, and 
infer that the bilayer spin correlations weaken with increase in hole doping. When the low-energy 
optic excitations are taken into account we measure the total integrated weight around (tt, tt), for 
energies below 120 meV, to agree with that expected from the insulator. As a qualitative guide, 
we compare spin-wave calculations for an ordered and a disordered stripe model and describe the 
inadequacy of this and other stripe models for the high-energy fluctuations. 

PACS numbers: 74.72.-h, 75.25.+Z, 75.40.Gb 



I. INTRODUCTION 

The discovery of superconductivity in the cuprates has 
resulted in a vast amount of interest in the spin fluctua- 
tions of these nearly two-dimensional materials. De- 
spite the intense research the cuprates have received, 
there is no generally accepted theory of the mechanism 
for the high temperature superconductivity. Due to the 
interplay between antiferromagnetism and superconduc- 
tivity the spin fluctuations are drastically altered upon 
doping the parent insulator and are also strongly corre- 
lated with superconductivity. Understanding the spin- 
fluctuations in the cuprates is an important step to un- 
derstanding superconductivity. 

We have recently completed an extensive experimental 
study of the temperature, momentum, and energy de- 
pendence of the spin fluctuations below 40 meV, includ- 
ing the resonance peak at ~ 33 meV, for the partially 
detwinned YBa2Cu306.5 ortho-II superconductornS, The 
good structural properties from the oxygen chain order 
removes the effect of impurity scattering of the quasipar- 
ticles and spin fluctuations due to disorder. We found the 
low-energy excitations to be gapless and to give rise to 
incommensurate peaks only displaced along the [100] di- 
rection from the (tt, tt) position. This contrasts with the 
symmetric cone of scattering which disperses from the 
(tt, tt) position observed in the parent insulating materi- 



als and illustrates the need for studies to be conducted 
on detwinned samples so that the a — b anisotropy can be 
unambiguously resolved. At ~ 33 meV we found a com- 
mensurate resonance peak which was strongly enhanced 
in the superconducting state. One of the surprising as- 
pects of the resonance peak is that, despite the low energy 
scale of superconductivity (Tc~59 K~ 5 meV), supercon- 
ductivity has a very strong effect on the high-energy ex- 
citations around the resonance. An important question 
is then how the excitations above the resonance energy 
are affected by superconductivity and hole doping. 

The nature of the high-energy spin fluctuations has 
previously been studied only in the parent insulating 
cuprate compounds. Using spallation neutrons Hay den et 
alA were able to map both the acoustic and optic modes 
in YBCOg.is and found that the results were well de- 
scribed by linear spin- wave theory. They found a large ~ 
75 meV energy gap for the optic fluctuations, an in-plane 
exchange constant of ~ 125 meV, and a bilayer coupling 
of 11 meV. The dispersion was found to extend up to 
~ 250 meV. These results have also been confirmed by 
reactor measurements by Reznik et al^ The entire spin- 
wave dispersion of the single layer compound, La2Cu04, 
has been studied in detail by Coldea et al.^ 

After this paper was submitted, there have been new 
developments because of letters to Nature by Hayden 
et air" and Tranquada et al& on the high-energy spin 
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fluctuations in YBCOe.e and LBCO. We have there- 
fore included discussion and comparison with the square- 
like pattern in these cuprates and will show that the 
spin fluctuations on ordered YBCOg.s are quite different 
and specifically, they are isotropic in wave vector about 
Q=(7r, ^). 

The high-energy excitations of a doped superconduct- 
ing sample were first measured by Bourges et alM' in a 
twinned crystal of YBCOg.s- They found a broadening 
in IQI of the correlated peak (at Q=(7r, tt)) with increas- 
ing energy. Based on this result it was suggested that 
the spin fluctuations disperse (at least at high-energies) 
in a manner similar to that of the parent insulating com- 
pound. The momentum dispersion of the magnetic ex- 
citations around the resonance was later studied in de- 
tail by Aral et ali^ on, again, a twinned YBCOg.v sam- 
ple. They found evidence for two modes which meet 
at the resonance energy, one of which opens downwards 
and gives the incommensurate scattering measured by 
other groupsfii and a new mode which opens upwards 
in energy. A similar result has also been suggested 
by a recent study on optimally doped YBCOg.gs by 
Reznik et alw^ Another, recent experiment, on twinned 
Lai.875Bao.i25Cu04 has found high-energy excitations 
which are consistent with the dispersion of two-leg lad- 
ders (when S(Q, lo) is averaged over domains)i& Because 
all of the samples studied have been twinned, it is not 
clear if the high-energy scattering results from incom- 
mensurate peaks displaced along a particular direction 
(as we have observed in our reactor low-energy study) or 
if the excitations form a symmetric cone which is simi- 
lar to that of the parent insulator (see the lower panel 
of Fig. 121 for an illustration). It is clearly important to 
compare the lineshape of the magnetic scattering above 
and below the resonance energy on a detwinned sample 
to resolve how the one-dimensional incommensurate scat- 
tering evolves to high-energies. Our present study on a 
largely detwinned sample enables the intrinsic nature of 
the magnetic scattering at high-energies to be determined 
free of scattering by structural disorder. 

The nature of the spin susceptibility above the reso- 
nance energy has many important implications for the 
theory of the cuprates. Recent effective band calcula- 
tions using the random-phase approximation predict the 
presence of spin excitations that have a velocity much 
higher than that of the parent insulator In contrast, 
linear spin-wave models based on a stripe ground state 
predict that the excitations at high energies are nearly 
symmetric around the (tt, tt) position and disperse almost 
linearly with energy^iii^iiSiii Quantum critical theories 
of two-dimensional Heisenberg antiferromagnets suggest 
that the resonance energy can be interpreted as a gap 
which softens as a critical point is approached. Such a 
theory would suggest that the excitations well above the 
resonance are qualitatively similar to those of the parent 
insulating compound^ Therefore, a detailed knowledge 
of the true spin spectrum (involving dispersion and line- 
shape) above the resonance energy is crucial in under- 
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FIG. 1: (Color online) The scattering diagram for a constant 
in-plane wave vector of (H,H) = (l/2,l/2) shows that as the 
energy transfer, E, varies, so does the c-axis component of 
wave vector L(E), thus enabling the acoustic (odd) and optic 
(even) symmetries of the spin response excitations to be de- 
termined by combining data. In the left panel the 150 meV 
incident neutron wave vector k, is brought in along the [0 1] 
direction. The primarily acoustic component of the response 
at E=77 meV (solid lines) is obtained near L=5 as well as 
the optic response at 49 meV near L=3. To obtain a difler- 
ent relation between E and L, so as to determine the ratio of 
acoustic to weight at different energies, the incident beam is 
brought in at variety of angles, ip, to the [0 1] crystal axis 
(right panel). 

standing the cuprates and in constraining theories. 

We present a high-energy neutron inelastic study on 
a well-ordered superconducting YBCO sample. We in- 
vestigate the momentum and energy dependence of the 
magnetic excitations above the resonance energy. This 
study complements our previous work on the low-energy 
spectrum using a reactor neutron sourcei^ The paper is 
divided into two sections and a discussion. The first sec- 
tion describes the experimental setup. The second sec- 
tion describes how we extracted the acoustic and optic 
weights using the known bilayer structure factors as well 
as the dispersion at high-energies. We discuss the inte- 
grated intensity in terms of the total moment sum rule. 
In the Appendix we describe linear spin-wave calcula- 
tions for an ordered and disordered stripe ground state 
and compare the results qualitatively to the experimental 
data. 



II. EXPERIMENTAL DETAILS 

The sample consisted of six orthorhombic crystals with 
total volume ^ 6 cm'^ aligned on a multi-crystal mount 
with a combined rocking curve width of ^ 1.5°. Details 
of the crystal growth, the detwinning by stress along the 
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FIG. 2: The upper panel shows the calculated acoustic struc- 
ture factor at the Q=(l/2, 1/2) position as a function of en- 
ergy transfer. It shows how one can choose a different frac- 
tion of acoustic and optic weight at at given energy transfer 
by changing the angle (t/i) between the incident beam direc- 
tion and the c-axis. The lower panel shows, schematically, 
the observed difference between scattering above and below 
the resonance energy. The scattering below the resonance en- 
ergy consists of two peaks displaced along the [100] direction 
while the high-energy scattering comprises a ring. The ar- 
rows illustrate how cuts were done to study the dispersion 
of the correlated peaks. The filled boxes (denoted A and B) 
represent the integration windows used to study the spectral 
weight and the lineshape as a function of energy transfer as 
discussed in the text. 



a axis, and the oxygen order have been published ear- 
lieri^i^ The sample has previously been measured to be 
partially detwinned with the majority domain occupying 
70 % of the total volume and with good oxygen order. 
Oxygen correlation lengths exceed ~ 100 A in the a and 
b directionsiSS The DC magnetization shows a supercon- 
ducting onset temperature of 59 K with a width of ~ 2.5 
K. The low temperature lattice constants were measured 
to be a=3.81 A, &=3.86 A, and c=11.67 A. 

Neutron scattering experiments were conducted by 
means of time-of-flight and triple-axis techniques. Time- 
of-flight measurements were made at the MAPS spec- 
trometer at the ISIS spallation source. The MAPS 
position-sensitive-detectors cover 16 m^ area divided 
into 36864 approximately square pixels. The position- 
sensitive-detectors are designed so that there is no di- 



rection that has significantly poorer Q resolution arising 
from the shape of the detectors. A monochromatic beam 
was produced by the use of a Fermi chopper spinning at 
a frequency of 400 Hz with the incident energy fixed at 
Ei=150 meV. The sample was mounted with the [001] 
and [110] directions in the horizontal plane and with the 
incident beam (ki) initially aligned along the [001] di- 
rection. Because the value of L, and hence the bilayer 
structure factor, varies with energy transfer (see FigQ] 
and the next section), this orientation gives only a lim- 
ited ability to measure the acoustic and optic response 
independently. In order to extract the acoustic and optic 
components for a range of energies (see Fig. 01, we varied 
the angle between the incident beam and the [001] axis, 
defined as ip, with the rotation axis along [110]. This pro- 
cedure was chosen over varying the incident energy (used 
elsewhere ''^'') to avoid changes in the resolution function. 
The full- width at half-maximum (FWHM) energy resolu- 
tion was measured from the incoherent elastic scattering 
to be ^ 8 meV. Using this method we have been able to 
study the spectral weight and the dispersion of both the 
acoustic and optic modes up to high-energies. 

To calibrate the spectrometer, a standard vanadium 
sample with a known mass was used. To obtain x"(q, uj) 
we have taken the cross section to be isotropic in spin 
with definitions of x" a-nd the total cross-section that 
we used previouslj*^ and that agree with those of other 
groupsiSi The definitions used here are identical to that 
used in Ref. ^ and in Ref. |3^ The definition of x" used 
here differs to that used in Ref. [s^ such that our values of 
x" are 1.5 times greater than that in Ref. js^- As a con- 
sistency check (and as discussed later) we have found that 
the total integral (/ d^q J dio) of the low-energy scatter- 
ing obtained from the vanadium calibration agrees well 
with that previously obtained'^ using an internal phonon 
calibration. 

Triple-axis measurements were conducted at the DU- 
ALSPEC spectrometer located at the C5 beam of the 
NRU reactor at Chalk River Laboratories. A focusing 
graphite (002) monochromator and a graphite (002) an- 
alyzer were used. A pyrolytic graphite filter in the scat- 
tered beam eliminated higher order reflections and the 
final energy was fixed at Ey=14.6 meV. The horizontal 
coUimation was set at [33' 29' S 51' 120']. The vertical 
coUimation was fixed at [80' 240' S 214' 430']. The six- 
crystal assembly was mounted in a closed-cycle refrigera- 
tor and aligned such that reflections of the form (HHL) 
lay in the scattering plane. To put these measurements 
on an absolute scale, a transverse acoustic phonon was 
measured around the (0 6) position. This was the same 
procedure used previouslyi^ 



III. RESULTS 

In magnetic systems that are two-dimensional (eg. 
La2Cu04), or one-dimensional (eg. CsNiCla), the non- 
dispersive nature of the excitations along one or two di- 
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rections makes the time-of-flight technique weU suited for 
mapping out the dispersion of the excitations. For such 
cases the non-dispersive momentum direction is typicaUy 
placed parallel to the incident wave-vector such that the 
scattering can be projected on the plane perpendicular 
to ki and the entire dispersion (in q^: and qy) can be 
mapped out. 

This technique has been successfully applied to the 
single-layer La2Cu04 insulator where the excitations are 
nearly two-dimensional and there is little dispersion along 
the [001] axis. The situation in YBa2Cu306+x is com- 
plicated by the presence of two strongly coupled Cu02 
layers in the unit cell. There is a weak antiferromag- 
netic coupling (of order 10 meV) between the two layers 
(separation d) in each unit cell. This means that the spin 
excitations can be categorized as odd (or acoustic) where 
the spin fluctuations are out of phase between the layers, 
or even (or optic) where they are in phasei^ Because the 
bilayer coupling is antiferromagnetic the optic spin waves 
in the insulator have a larger energy than the acoustic. 

The acoustic mode intensity varies along the [001] as 
sm^{Qzd/2), where d is the distance between the two 
layers in the unit cell, while the optic mode strength is 
proportional to cos'^{Qzd/2).^ Therefore, for a fixed ex- 
perimental configuration different L values, and hence 
different acoustic and optic weights, will be sampled be- 
cause L varies with energy transfer. The reason for L 
varying with energy is schematically outlined in Fig. ^ 
which shows the scattering triangle for two different en- 
ergy transfers, and hence two different final energies (note 
that E, is fixed in this experiment). For a particular 
energy transfer corresponding to a final energy Ej, the 
constraint of having the momentum transfer of the form 
Q=(l/2, 1/2, L) fixes the value of L. Therefore, by chang- 
ing the energy transfer, but keeping Q to have the form 
Q=(l/2, 1/2, L), will result in L varying with energy 
transfer. 

The effect of this is further illustrated in the upper 
panel of Fig. 13 which shows the acoustic structure fac- 
tor at Q=(l/2, 1/2) as a function of energy transfer for 
■0=0° and 20°. The acoustic structure factor varies along 
[001] as sin'^{Qzd/2), but since the value of Qz depends 
on the energy transfer for a fixed {Qx,Qy)={l/2, 1/2), 
different acoustic and optic weights will be sampled as 
a function of energy. Therefore, for tp^O (see Fig. 
to sample purely acoustic scattering an energy transfer 
of about 30 meV or 75 meV is required. Fig. [21 also il- 
lustrates that varying moves the peak of the acoustic 
structure factor to different energy transfers. In other 
studies, this problem has been overcome by tuning the 
incident energy so that the energy position of the maxi- 
mum in the structure factors would shift to the required 
position. Such a procedure inevitably alters the reso- 
lution function and therefore the relative weight of the 
acoustic and optic modes. As a result, the dispersion of 
the acoustic and optic modes cannot be reliably deter- 
mined without the aid of resolution corrections. 

To avoid this, we kept the incident energy fixed and 



varied the angle between the [001] axis and k.^ (il- 
lustrated in Fig. P). We have collected data for ip=0, 
20°, and 30°. In each configuration we were able to ex- 
tract data in up to 12 different Brillouin zones (including 
the (±l/2,±l/2), (±3/2,±l/2) and permutations) each 
of which occur for different values of L. The intensity of 
the magnetic scattering for each zone, sample orientation 
("0), and energy bin was then obtained by taking a cut 
(along the [1, 0] and [0, 1] directions) through the scat- 
tering at a constant energy transfer. This process relied 
on the MSLICE program,^^ 

In this section we discuss the experimental results in 
terms of four different subsections. In the first section we 
discuss the lineshape of the high-energy acoustic scatter- 
ing and its symmetry in the (H, K) plane. This is based 
on two-dimensional slices integrating ±7.5 meV in en- 
ergy. The next section we study the momentum depen- 
dence of the scattering and its dispersion in detail. This 
involves an analysis where we integrated over fine slices 
in momentum, (see Box B in the lower panel of Fig. O 
and therefore do not smear out any lineshape changes 
as a function of energy transfer. In the third section 
we describe the relative optic and acoustic weights as a 
function of energy transfer using an analysis in which we 
integrate over one momentum direction. This is schemat- 
ically represented by Box A in the lower panel of Fig. |2 
This removes any detailed lineshape information and al- 
lows all energies to be treated equally. In the final section 
we discuss the total integrated spectral weight and com- 
pare the results to other systems as well as the insulating 
compound (in particular YBCOg.is). 



A. Spin Excitation Ring at High-energies 

Examples of two-dimensional constant energy maps 
taken at ^—0 are shown in Fig. 13 The two-dimensional 
slices shown in Fig. Ohavc been averaged over the four (± 
1/2, ± 1/2) zones sampled by the MAPS detector array. 
At the resonance energy of 33 nieV the response is local- 
ized at (tt, tt). Above the resonance at 78 meV a ring is 
observed showing the intersection of the constant energy 
surface with a two-dimensional cone of spin waves cen- 
tered on (tt, tt) (this point is discussed more thoroughly in 
the next section). The intensity around the ring is sym- 
metric in the (II,K) plane as demonstrated in Fig. ^ Fig. 
^ plots one-dimensional cuts through the (tt, tt) position 
at E=77 meV along the [1, 0], [0, 1], [1, 1], and [1, T] di- 
rections. For comparison purposes a sloping background 
has been subtracted from each cut. The solid lines are 
Gaussians symmetrically displaced from the (tt, tt) posi- 
tion. A comparison between the solid lines and the data 
shows that the scattering is consistent with an isotropic 
ring. Both the radius on wave vector and the intensity 
are independent of the direction in the a-b plane. 

Maps at all energies higher than the resonance are 
consistent with almost isotropic spin waves centered at 
(tt, tt). Although an individual constant-energy contour, 
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FIG. 3: (Color online) Smoothed two-dimensional slices 
through the correlated response for ip^O at 32.5 meV and 
77.5 meV energy transfers and integrated it 7.5 meV along 
the energy axis. At 77.5 meV the ring of scattering shows 
the intersection with a constant energy surface of the cone of 
spin wave dispersion emanating from the (1/2, 1/2) position. 
Within statistics the velocity is isotropic. The intensity rep- 
resented by false color is in arbitrary units. Unsmoothed cuts 
of the data for various energy transfers are presented in Fig. 
Hand Fig. M 



such as that at 77 mcV in Fig. |21 may show more counts 
at particular directions around the ring, our results are 
statistically consistent with a spin- wave dispersion whose 
velocity and spectral weight are isotropic in the (H, K) 
plane. The line shape illustrated in the color map in 
Fig. O may also be compared with the diamond shape 
suggested by various calculations involving ladders 
For the YBCOe.s system the scattered intensity clearly 
does not correspond to a square pattern nor to peaks 
displaced along [110] (see Fig. There is no evidence 
at energies greater than 33 meV for the one-dimensional 
incommensurate modulation that exists below the reso- 
nance energy^ 

Our observation that the scattering is symmetric is im- 
portant in the context of a recent model proposed by 
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FIG. 4: One-dimensional cuts through the correlated ring 
at 77 meV energy transfer (displayed in Fig. O along the 
[10], [01], [11] and [11] directions. The data was integrated 
± 0.05 r.l.u. perpendicular to the cut direction and a sloping 
background was subtracted from each scan for comparison 
purposes. The solid lines are Gaussians displaced equally from 
the (tt, tt) position. 



Tranquada et al. which suggests that the high-energy 
scattering in Lai.875Bao.i25Cu04 (LBCO) is similar to 
the excitations from ladders At high energies compa- 
rable to ours, the actual data in the color plots of Ref. 
1^ seem to the eye to show less of the strong diamond- 
shaped anisotropy of the ladder model and are not dis- 
similar to the symmetric excitations we observe. We are 
able to resolve the excitation in wave vector because of 
the lower velocity in YBCOe.s. 

Recent calculations by Vojta and Ulbricht have found 
that a bond-centered stripe phase has an excitation spec- 
trum very similar to that observed in LBC0i2^ For a 
ladder ground (or bond-centered stripe) state to give a 
symmetric excitation spectrum, equal contributions from 
horizontal and vertical ladders (or stripes) must be con- 
sidered, as is present in a twinned superconductor. Our 
results on partially detwinned YBCOg.s show that a de- 
twinned structure does not give rise to a large anisotropy 
such as that expected for ladders. This anisotropy should 
be particularly evident when comparing cuts along [11] 
to those along [10] or [01]. Ref. H modelled the ladders 
for a twinned system (ladders along a and b) and so may 
give stronger scattering along [11] due to the overlap of 
the a and b domains. It may be that the stripes (or 
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ladders) are purely dynamic and therefore higher energy 
fluctuations do not sense the small orthorhombic nature 
of the crystal structure, even though the low energy fluc- 
tuations do. We note that in the energy range studied in 
our experiment the dispersion relation for a ladder (like 
that observed in LBCO) is similar to the linear spin- 
waves observed in the insulator. 

We will show later that the high-energy excitations 
have a similar dispersion and spectral weight to those of 
the parent insulator YBCO6.15. Therefore, incommen- 
surate fluctuations that are displaced and are stronger 
along a particular direction occur only at low-energies 
indicating that the stripe texture barely influences the 
high-energy spin excitations. This suggests that the high- 
energy excitations exhibit similar isotropy to those of an 
undoped insulator. 

Our results differ from those reported in a recent study 
on twinned YBCOg.e by Hayden et al7 who measured, 
not a ring of scattering, but four peaks displaced along 
the [110] directions. For the YBCOg.s ortho-II detwinned 
system the scattered intensity clearly does not corre- 
spond to a square pattern nor to peaks along [110] (see 
Fig. ^ . The different behavior may arise from the larger 
x=6.6 doping, for the low-energy incommensurate wave 
vector {5) in Ref. is at least twice as large as that in 
Ortho-II YBCO6.5. Since 5 scales approximately with 
doping, the YBCOe.e system of Ref. [3 must have a sig- 
nificantly larger effective doping in the Cu02 plane than 
our Ortho-II YBCOg.s. We note that modestly larger Tc 
of 63 K'' is not a direct measure of doping since it is, in 
general, lowered by increased structural disorder. 



B. Dispersion of the Optic and Acoustic Modes 

Despite the fact that the low-energy excitations be- 
low the resonance energy show one-dimensional incom- 
mensurate peaks slightly displaced along the [100] direc- 
tion, the high energy excitations are close to being sym- 
metric around the (tt, it) position. This is illustrated 
in Fig. |31 which shows two-dimensional constant energy 
slices through the correlated peaks at T = 6 K. In par- 
ticular, there is no preferred weight of the spin response 
at incommensurate wave vectors along a* relative to b* . 
Any modulation around the spin-wave cone lies within 
statistical error as demonstrated by the one-dimensional 
cuts displayed in Fig. ^ Thus, the high-energy scattering 
forms a symmetric cone, whose momentum increases with 
increasing energy as expected for the linear, constant- 
velocity part of the spin-wave dispersion. We have fur- 
ther observed that there is no qualitative change in the 
momentum dependence of the high-energy scattering in 
the normal state at T=85 K. 

To investigate how the lineshape evolves as a function 
of energy transfer, one-dimensional cuts along the [100] 
axes at each energy transfer were fitted to two Gaussians 
displaced equally from the (tt, tt) position to obtain the 
dispersion for the acoustic and optic modes above the res- 
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FIG. 5; The dispersion of the acoustic and optic modes with 
respect to the (tt, tt) position at T= 6 K is plotted. The filled 
circles represent the peak positions of a two Gaussian fit to 
cuts along the [100] and [010] directions above the resonance 
energy and from cuts along only the [100] direction below the 
resonance energy. The open circles represent the positions of 
the incommensurate peaks found in experiments conducted 
at Ghalk River. The open squares represent the peak posi- 
tion that result from a polar average around the {n,n) posi- 
tion. The solid lines schematically represent the dispersion of 
the insulating compound as measured by Hayden et al. The 
dashed lines are fits of the high-energy dispersion to linear 
spin- wave theory for Tluj > 40 meV as described in the text. 
The horizontal bars show the q width observed for three optic 
scans at constant energy. 



onance energy. As a consistency check, we have also per- 
formed the same analysis above the resonance but with 
cuts along [010]. Identical results were obtained. The 
results for the overall momentum dependence of the spin 
fluctuations is illustrated in Fig. |S1 Unlike the analysis 
described in the next section, the one-dimensional cuts 
to investigate the dispersion involve integrating over only 
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FIG. 6: Constant energy cuts through the correlated peak 
at the {n,n) positions are plotted. An energy integration of 
±7.5 meV was used and the data are integrated ±0.05 r.l.u. 
along the [010] direction. The solid lines are guides to the 
eye. By conducting similar constant energy cuts to those 
displayed here the dispersion curve was obtained up to ~ 100 
meV energy transfer. 



a small range perpendicular to the cut direction. This 
procedure is schematically illustrated in Fig. [21 by the 
integration window represented by Box B. Examples of 
such one-dimensional cuts are shown in Fig. for the 
acoustic mode and Fig. [7|for the optic mode. 

For energies below the resonance we have confirmed 
that the stronger incommensurate scattering, along [100], 
previously obtained with a triple-axis spectrometer is 
consistent with the MAPS data. The low-temperature 
results previously obtained using a triple-axis spectrom- 
eter (and described elsewhere) are plotted as open cir- 
cles in the upper panel of Fig. [S] The result at T=6 
K is plotted in Fig. [S] Here we have imposed a cut- 
off for each mode such that acoustic is defined as any 
scan with sm^{Qzd/2) > 0.8 and likewise for the optic 
mode cos^{Qzd/2) > 0.8. Based on this strict criterion 
we expect the measured dispersions to reflect the acous- 
tic and optic dispersion accurately. The data presented 
in Fig. 13 were obtained for an energy integration width 
of AE=10 meV and an integration along [010] of ± 0.05 
r.l.u. We have also verified that the overall results re- 
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FIG. 7: Constant energy cuts through the correlated peak 
at the (7r,7r) positions are plotted. An energy integration of 
±7.5 meV was used and the data are integrated ±0.05 r.l.u. 
along the [010] direction. The solid lines are guides to the 
eye. By conducting similar constant energy cuts to those 
displayed here the dispersion curve was obtained up to ~ 100 
meV energy transfer. 



main unchanged for larger energy bins. The horizontal 
errorbars on each point are the errors in the q positions, 
for two branches except at the resonance energy where 
only one peak is resolved. The open circles in the up- 
per panel of Fig. [S] are locations of the incommensurate 
peaks obtained from our previous reactor experiment on 
the acoustic mode. The reactor results are discussed in 
detail in Ref. 0, and hence will not be presented in detail 
here. Despite the fact that the energy resolution used in 
the current experiment is much broader than that used 
in our reactor work the two results agree very well. 

The open square symbols in Fig. El represent a polar 
average around the (vr, tt) position obtained for scans 
done with ijj=0, i.e. with the incident beam parallel to 
the [001] direction. For all energy transfers in the range 
huj=A5 meV to 105 meV, in steps of 5 meV, the polar 
average was evaluated with respect to the (tt, tt) position. 
Points were binned in 0.017 rlu steps from (tt, tt). For 
the points in each bin the mean distance from (tt, tt) and 
the mean intensity were evaluated. A Gaussian fit to 
the resulting polar average was used to find the position 
of the maximum intensity for that energy transfer. The 
polar average agrees well with the one-dimensional cuts 
along [100] and [010] as we expect for an isotropic velocity 
for the high-energy excitations. 
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FIG. 8: A map of the acoustic intensity (color online) as 
extracted from Gaussian fits to constant energy scans. The 
intensity has been plotted on a logarithmic scale and is in 
arbitrary units normalized to unity. From this plot, it can be 
seen that an energy scan at constant Q, at energies exceed- 
ing the resonance energy, will not display a peak or evidence 
of propagating modes. There is nonetheless evidence of two 
ridges that move to larger momentum as energy increases and 
that give a qualitative indication of a characteristic velocity 
(dotted line). 



The momentum dependence of the spin fluctuations as 
a function of energy transfer qualitatively differs from the 
scattering measured in early work on Lai.86Sro.i4Cu04 
which appears to show no appreciable broadening in mo- 
mentum up to energies of ~ 150 meV^ Our results 
seem to suggest that the resonance may be interpreted 
as the top energy scale of the stripe-like incommensurate 
features at low-energies. This is confirmed by a recent 
experiment on Lai.84Sro.i6Cu04 by Christensen et al.-'^ 
They found that the low-energy incommensurate peaks 
do meet at the (tt, tt) position at ^ 40 meV, a simi- 
lar energy scale to that found in YBCO. As suggested 
in Ref. |23j the similarity between this and the 41 meV 
resonance measured in optimally doped YBCO may sug- 
gest that the resonance is not directly associated with Tc- 
The strong similarity between the low-energy response in 
LSCO and YBCO may imply that the spin spectrum is 
common to all cuprates. The low-energy responses mea- 
sured in ortho-II YBCOg.s and Lai.84Sro.i6Cu04 are also 



qualitatively similar to the spin spectrum found in the 
nickelates, a well characterized system known to exhibit 
charge and spin stripes iSSiSS 

The dispersion measured in YBCOg.s contrasts with 
that of the insulating compound. The solid lines in 
Fig. 13 represent the dispersion measured in insulating 
YBCOe.ig^. For the acoustic modes the lower dashed 
lines represent a fit to linear spin-waves {huj oc cq) to 
only the highest energy excitations {huj > 40 meV). For 
the optic modes, the dashed line is the result of a fit to 
linear spin-wave theory for bilayers. We obtain a slope of 
he ^ 400 meV • A for the high-energy acoustic scatter- 
ing in our well-ordered ortho-II crystal. The spin-wave 
velocity observed in superconducting YBCOe.s is consid- 
erably reduced from the insulator value of he ^ 650 meV 
• All It is also considerably less than the value of he ~ 
600 meV • A measured recently for LBCOi^ 

The acoustic dispersion around the resonance energy 
is very similar to that to the a;-like dispersion observed 
in Lai.875Bao.i25Cu04)^ LSCO)Si and optimally doped 
YBCOji^ This suggests that the a;-like dispersion is com- 
mon to all cuprates. Despite the fact that the disper- 
sion is similar, our results in detwinned and oxygen or- 
dered YBCOe.s show the low-energy excitations to be 
one-dimensional with incommensurate peaks displaced 
only along the [100] directions and the high-energy ex- 
citations to be essentially symmetric around the (tt, tt) 
position. We note that a similar value to that obtained 
here for the spin-wave velocity was found by Bourges 
et alA for disordered superconducting YBCOg.s {T^c— 
52 K). Therefore, the high energy excitations are not 
much influenced by structural disorder, whereas disorder 
broadens the intrinsic sharpness of the resonance consid- 
erably.^ This emphasizes our previous assertion that the 
high energy excitations may not be that sensitive to the 
effects of twinning and the orthorhombic anisotropy. 

We note that an important difference between the work 
presented here and that of previous studies (such as Ref. 
|0) is that by combining reactor and spallation data we 
have been able to resolve the dispersion around the reso- 
nance and contrast the high-energy spin excitations with 
the incommensurate rods at low-energies. Previously, it 
was not clear if the dispersion at high-energies extrapo- 
lated smoothly (as interpreted in Fig. 3 in Ref. 0) to 
huj=Q or not. Our results clearly show that the disper- 
sion at high-energies is distinctly separate from the low- 
energy incommensurate fluctuations. The separation be- 
tween the two regions is defined by the resonance energy. 
The physical mechanism separating these energy regions 
is not currently clear, however, the high-energy excita- 
tions are of a different character than those observed at 
low-energies below the resonance. 

To estimate the in-plane (J||) and out-of-plane (J±) 
couplings, we have fitted the acoustic dispersion from 
the resonance to high energies (ignoring the low-energy 
incommensurate scattering) to 

EacousUc{q)^{^L + <qfY'\ (1) 
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where e{q) = ficq (at low-energies) is the acoustic dis- 
persion predicted by Hnear spin-wave theorji^iS^ and an 
arbitrary gap A^c ~ 33 meV has been set close to the res- 
onance energy. For the optic mode, we used the gapped 
form predicted by linear spin-wave theory and found Aop 
~ 50 meV from the highest -energy optic modes (lower 
panel of Fig. (S)). Based on this analysis we have obtained 
J||'^70 meV and 10 meV. The derived velocity ~ 365 
meV • A, using the gapped form in Eq. ^ is close to that 
obtained by extrapolating the high-energy spin-waves to 

The use of the gapped form for the acoustic disper- 
sion is physically justified by the fact that the low-energy 
incommensurate branches meet at the resonance energy 
therefore suggestive that the low-energy incommensurate 
scattering and the high-energy dynamics can be treated 
as different branches. Also the gapped dispersion results 
naturally if we consider the ground state to be approxi- 
mated by dynamic stripes or ladders.'^" 

It is well known that for the 5"= 1/2 square lattice, 
renormalization factors are important. To estimate the 
true exchange constants the renormalization factor Z^, 
has previously been used. This has been estimated from 
a 1/S' expansion by Igarashi to be ^ 1.2m^ Even though 
it is only a 20 % effect, the effect of the renormaliza- 
tion factor means that our measured value of 70 meV 
translates into ^ 60 meV. The value of J|| of 60 meV is 
considerably less than the value of 125 meV found in the 
insulating compounds. We conclude that the effect of 
doping has been to decrease the spin exchange, J||, and 
concomitantly, the velocity by 50 %. 

For energies below ~ 40 meV we were not able to sepa- 
rate reliably the two peaks in q associated with the optic 
mode. The dashed bars in Fig. ISjrepresent the full width 
at half maximum AqopUc of the optic spin response at ~ 
30 meV. We believe the incomplete spin correlations (see 
next section) may be responsible for the optic tail be- 
low ~ 40 meV, and so the data at low energies should 
not be fitted to our model dispersion. If the optic tail is 
the result of dynamic stripes then the strong similarity 
between the high-energy dynamics from both the disper- 
sion and the integrated intensity (see next section) sug- 
gests that the stripes are weakly pinned and therefore the 
high-energy scattering can be treated differently than the 
scattering around, or below, the resonance energy. Our 
suggestion of the presence of weakly pinned stripes is con- 
sistent with recent resonant x-ray scattering data illus- 
trating a weak hole modulation in the Cu02 planes, pos- 
sibly associated with chain ordering. '^^ If the stripes are 
weakly pinned, then the spin fluctuations at low-energies 
would be expected to have very different properties from 
those at high-energy. This provides some support for our 
analysis in which we extract effective exchange constants 
from only the high-energy data. 

A prediction of the bilayer stripe modeU^ is that the 
optic scattering near the resonance energy will be in- 
commensurate. It is difficult to discuss the possibility 
of incommensurability in the optic channel in YBCOg.s 



because the optic scattering is relatively weak and be- 
cause we knowS^ that the incommensurate wave vector is 
small, 6 ~ 0.06 r.l.u.. As shown in Fig. [S] the width 
of the optic scattering below ~ 40 meV does not change 
measurably when analyzed as a single peak. However the 
scan at 32.5 meV in Fig. [7|a gives a slight hint that the 
low-energy optic scattering may also be incommensurate 
with a S similar to that of the acoustic channel. Further 
studies at higher doping where S is larger may be useful 
in resolving this issue. 

Although constant energy cuts show well defined 
peaks, we find that constant-Q cuts up to ~ 100 meV 
do not show any clear sign of propagating modes. Fig. 
IHl plots the fitted intensity and widths as a function of 
energy transfer based on constant energy cuts through 
the data from 20 to 90 meV. Based on this figure, it can 
be seen that constant-Q scans do not show evidence for 
propogating modes (ie. a peak in a constant-Q scan) 
above the resonance energy. This may suggest that the 
spin fluctuations are heavily damped even at high ener- 
gies but to confirm this our results need to be extended 
to much higher energies due to the large velocity of the 
spin excitations in the energy range studied here. It is 
interesting to note that despite the fact we observe two 
well defined peaks in the acoustic channel (Fig. it is 
difficult to resolve the two optic modes even at ~ 93 meV 
(Fig[7)). This may indicate that the optic channel is more 
heavily damped both in momentum and energy, than the 
acoustic channel. 



C. Extraction of Acoustic and Optic Intensities 

The previous sections have illustrated that the line- 
shape of the magnetic scattering varies dramatically 
with energy transfer. When comparing relative spectral 
weights between the acoustic and optic modes, it is desir- 
able to be able to compare each energy equally without 
need for detailed lineshape corrections. In this section 
we discuss our analysis of the acoustic and optic weights 
by integrating over one dimension in q. 

To analyze the magnetic intensity into its optic and 
acoustic components as a function of energy transfer we 
have assumed a symmetric two-dimensional Gaussian line 
shape (along [1, 0] and [0, 1]) and integrated over a broad 
range in |q|. The range is chosen wide enough to collect 
intensity from all the incommensurate features. This pro- 
cedure is illustrated in Fig. [21 by the integration window 
represented by Box A, which shows a broad integration 
along the [0, 1] direction perpendicular to the cut direc- 
tion. By conducting a one-dimensional cut along [1, 0], 
with the intensity integrated along [0, 1], we obtained a 
single peak that includes all of the correlated scattering. 
The single correlated peak was fitted to a Gaussian plus a 
sloping background. The width of the integral along the 
[0, 1] was taken to be ± 2a of the fitted Gaussian along 
[1, 0]. We repeated the analysis by doing a cut along [0, 
1] and integrating by ± 2cr along [1, 0] and obtained con- 
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FIG. 9: Scans through the acoustic correlated peak for several 
energy transfers. The scans have been integrated perpendic- 
ular to the scan direction over a broad range as discussed in 
the text in order to include all of the correlated scattering. 
The energy integration width for each scan is ± 2.5 meV. 



sistent results. Examples of this analysis are displayed 
in Figs. El and ^1 for scans of predominantly acoustic 
and optic weight respectively. The integrated intensity 
profiles at all energies are reasonably described by a sin- 
gle Gaussian, and more importantly, even when some of 
the underlying two-peaked structure remains visible, the 
Gaussian does a good job of estimating the integrated in- 
tensity needed as input to the ensuing the acoustic-optic 
analysis. 

The large integration width perpendicular to the scan 
direction has smeared out the ring of scattering (as shown 
in Fig. ISJ and allows the correlated intensity to be well 
described by a single Gaussian from which an amplitude 
can be extracted. The use of a one-dimensional integral 
to study the intensity as a function of energy transfer has 
the advantage of removing any changes in the line shape 
with energy transfer (as discussed in the next section). 
With this analysis we can treat all energy transfers on an 
equal footing in terms of a single symmetric Gaussian, 
and study the acoustic and optic components without 
need for incorporating detailed changes in the line shape. 
We discuss the two-dimensional integral as a function of 
energy transfer later. 

To extract the acoustic and optic weights at each en- 
ergy transfer, the amplitude from the Gaussian fit was 
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FIG. 10: Scans through the optic correlated peak for several 
energy transfers. The scans have been integrated perpendic- 
ular to the scan direction over a broad range as discussed in 
the text. The energy integration width for each scan is ± 2.5 
meV. 



corrected for the integration perpendicular to the scan 
direction and the anisotropic magnetic form factor. We 
emphasize that this analysis relies on the assumption of a 
symmetric Gaussian centered at the (tt, tt) position. The 
measured peak intensity of a cut along a particular di- 
rection through the peak of a two-dimensional Gaussian 
will depend on the integration window perpendicular to 
the cut direction. However, once the width of the Gaus- 
sian is known (defined by cr), the correction for the finite 
integration perpendicular can be computed to obtain the 
true peak intensity of the Gaussian. To check this analy- 
sis, several integration widths were chosen and consistent 
results were obtained. By plotting the corrected inten- 
sity as a function of the different L values we obtain the 
modulation data shown in Fig. ^] The depth of the 
modulation reflects the acoustic-optic ratio. The known 
structure factors, 



x"(Q.,c.) = x'L. sin'{Q,d/2) + Xop. cos2(Q,d/2), (2) 

then enable the acoustic [xlc ) ^-^id optic (x"p ) suscepti- 
bilities to be obtained from a least-squares fit to the data 
shown in Fig. In the fits we have fixed the parameter 
d to be the known bilayer spacing22. as successfully done 
at 33 meV in our earlier reactor experiment. We empha- 
size that the separation of spectral weight into acoustic 
and optic components is solely based on the fact that 
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FIG. 11: The intensity as a function of L is plotted at three 
different energies. The solid line is the result of a fit to the 
sum of the acoustic and optic structure factors described in 
the text. Based on similar fits, optic and acoustic intensities 
were extracted as a function of energy transfer as displayed in 
Fig. 1121 For the fitted curves, d/c was fixed at 0.3. The dotted 
vertical line shows where we scanned using a triple-axis to 
confirm the presence of optic spectral weight at low-energies. 



300 



> 

CD 

:x 200 

s 
o 

100 



YBCO 
T=6 K • 

r \ 


05 

(Chalk River) 


r ' X 




/ ' T 

m ' A 


o Optic 




ooo 




40 60 80 
E (meV) 



120 



FIG. 12: The intensities of the acoustic and optic modes are 
plotted as a function of energy transfer at T=6 K in the up- 
per panel. The solid bar represents the FWHM resolution of 
MAPS measured at Tiijj=Q. The acoustic and optic intensities 
are roughly equal at high-energy transfers. The optic mode 
can be seen to extend to very low-energies compared with the 
parent antiferromagnetic insulator. The solid lines are guides 
to the eye. The dashed line is the low-temperature lineshape 
of the resonance obtained using a triple-axis spectrometer and 
scaled down from ~ 1400 /ig / e V to 300 y?g /eV to account for 
the poorer energy resolution of the MAPS data. The lower 
panel shows the energy dependence of the average FWHM of 
one-dimensional cuts along H and K for which the transverse 
momentum integration (along K and H) is set large enough 
to cover the diameter of the spin response. 



there are two coupled layers per unit cell, and does not 
depend on correlations or damping of the two modes. 

Fig. displays the L dependence of the scattering at 
three energy transfers, because the L-dependence of the 
data extends to relatively low L, we have checked that 
the optic /acoustic ratio is insensitive to the precise value 
of d/c.The solid line is the result of a fit to Eq. 13 At low 
energies the scattering can be seen to be predominately 
acoustic while at high energies (above ~ 60 meV) the 
acoustic and optic modes contribute nearly equally at all 
L values. The nearly equal contribution for optic and 
acoustic modes is consistent with the results of Pailhes 
et al. at 54 meV in YBCOe.ss'^ The fact that we clearly 
see scattering out to L~ 7 suggests that the scattering 
is predominately from Cu^"*" spins. Some other forms of 
magnetism, like the orbital currents described in Ref. Is^ . 
would only contribute significantly at small values of L 



due to their steep form factor decrease. 

The results for the absolute intensities of the acoustic 
and optic modes are summarized in Fig. 1121 The ab- 
solute intensity as a function of energy transfer is qual- 
itatively consistent with data published previously for a 
similar oxygen concentration^ The strong asymmetry 
with a sharp cut-off at 35 meV that was observed in 
our reactor data (broken line in Fig. I12|) is observed less 
clearly because of the broader energy resolution of MAPS 
near the resonance. The lower panel of Fig. Elshows the 
FWHM, averaged over direction and zone, of a fit of a 
single Gaussian to the one-dimensional cuts at constant 
energy described previously. This figure is a qualitative 
picture of the dispersion already presented earlier in this 
paper. The FWHM gives the overall Q-width of the in- 
commensurate scattering at low-energy, the width of the 
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FIG. 13; Scans of the optic spin response at L=6.7 with the 
DUALSPEC triple-axis spectrometer at NRU reactor, Chalk 
River. It can be seen that optic intensity is still present at 
relatively low energy transfers. The solid lines are Gaussian 
fits to the data. The solid and open circles represent data 
taken at 9 K and 200 K, respectively. The 200 K has had 
200 counts subtracted to allow a direct comparison to the 9 
K data. 



commensurate resonance at the 33 meV resonance en- 
ergy, and at larger energies gives a qualitative measure 
of how the cone of spin excitations opens with increas- 
ing energy (discussed in detail in the previous section). 
The integrated data presented here (peak and FWHM) 
are used in the next section to compute the total integral 
overall momentum and energy. 

The acoustic fluctuations exhibit the strong resonance 
at 33 meV seen in our reactor experiments and also ex- 
hibits a weaker second peak at around 55 meV. This 
weak second peak may be consistent with the shoulder 
previously observed by Bourges et alA above the reso- 
nance energy. We note that a recent phonon study by 
Pintschovius et al.^^ has found intense and nearly disper- 
sionless branches which intersect with the zone boundary 
in this energy range. Therefore, we cannot rule out that 
any extra scattering may be the result of a significant 
phonon contribution in the range 50-60 meV. We note 
that we do not observe any new resonant feature above 
the 33 meV resonance, apart from the weak 55 meV fea- 
ture. Such a resonance has been suggested in theory to 
result from an antibound state of two holes, the tt par- 
ticlei^ If this new resonance is of a similar strength to 



that of the 33 meV resonance, our results put a lower 
bound of ~ 100 meV for the energy of any such resonant 
feature. Lee and Nagaosa^i have predicted that orbital 
current fluctuations cause additional scattering above the 
main resonance. Due to the form factor associated with 
orbital currents such scattering should only be apparent 
for small values of L, not accessible in our experimental 
configuration. It is currently not clear if orbital currents 
are contributing to the scattering we observe around 50- 
60 meV. 

A particularly interesting result displayed in both Figs. 
111! 1121 and ^1 is the presence of optical intensity at rel- 
atively low-energies. Our measurements show that if an 
optic gap exists, its energy is less than ~ 20-25 meV. Sur- 
prisingly, this is much lower than the optic energy gap of 
~ 75 meV measured in the insulating YBa2Cu306.i5>^ 

We have confirmed that the optical response extends 
to low energies by direct measurement with triple-axis 
spectrometer where both L and huj can be tuned to the 
desired values. Scans at the minimum of the acoustic bi- 
layer structure factor, L=6.7, are shown in Fig. Elat en- 
ergy transfers of 33.1 and 37.2 meV and at temperatures 
of 9 K and 200 K. The fact that the scattering decreases 
when we heat to 200 K supports the interpretation that 
the scattering is magnetic and not from phonons. The 
scans confirm our findings from MAPS that a correlated 
peak of optic symmetry exists at low energies at low tem- 
peratures. These are the same optic fluctuations that ex- 
tend to 25 meV in Fig. El A further consistency check 
is that at 33 meV the ratio of optic to acoustic weight ob- 
tained in our reactor experiment is consistent with that 
found in the MAPS data. The absolute intensity from the 
reactor data at 33.1 meV and 37.2 meV was estimated 
to be 142 ±16 ^J.%/eV and 164 ± 15 ^J.%/eV respectively 
based on a phonon calibration. This is compared with 
a peak value of ~ 1400 fj,%/eV previously measured for 
the acoustic resonance at 33 meV. The factor of about 
10 reduction from the acoustic resonance to the optic 
weight is similar to that measured at MAPS. The cor- 
related peak observed in our reactor experiment cannot 
arise from the feed through of the acoustic resonance at 
both energies, for we have previously demonstrated that 
the acoustic weight decreases substantially from 33 to 37 
meV. However, both our MAPS and DUALSPEC data 
shows the optic weight to be comparable. Also, previ- 
ous measurements of the resonance along [001]'^^ suggest 
that at L=6.7, there is essentially no acoustic contribu- 
tion to the scattering (or at least an attenuation factor 
much larger than 10). 

The Chalk River and MAPS data show that the optic 
fluctuations do not display any resonant behavior around 
33 meV unlike the strong acoustic resonance. A similar 
non-resonant optic response was obtained by Fong et al. 
in underdoped YBCOg+s who observed that the onset of 
the optic mode energy (~30 meV in YBCOe.?) declined 
with doping. '^^ Our results indicate the gradual build- 
up of non-resonant optic intensity around the resonance 
energy as has been observed for nearly optimally doped 
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The presence of optic weight at low-energies is surpris- 
ing when viewed in the hght of Hnear spin-wave theory 
which predicts a large optic gap of order ^ 70 meV for 
the antiferromagnetically ordered insulator. It is also a 
departure from models for the superconductor, such as 
the renormalized mean field theory of Brinckmann and 
Lee^, which predicts an optic gap equal to or larger than 
the acoustic resonance energy. In the previous section, we 
found that the measured high-energy dispersion requires 
a finite bilayer coupling of order 10 meV, which should 
have served to preserve the antiparallcl spin correlations 
between layers. The presence of optic spectral weight 
below the acoustic resonance, however, leads to the idea 
that there is a reduction in the correlation of the spins 
between the two layers. In the limit of strong disorder the 
correlation length along [001] is less than the bilayer spac- 
ing and there would be no L-dependence of the scatter- 
ing, analogous to the monolayer system. For weak disor- 
der the correlation length would exceed the bilayer spac- 
ing and the scattering would depend strongly on L. For 
YBCO an intermediate situation prevails. We speculate 
that the presence of optic spectral weight at low-energies 
signifies that the spin corelations in each layer are becom- 
ing increasingly weakened with increasing doping, hence 
reducing the modulation in their acoustic/optic spin am- 
plitudes. Strong decorrelation of the spins along [001] 
with increasing hole doping has been observed in both the 
La2-2:Sr2,Cu04 and La2Cu04+y systems in the elastic 
channel42iii The increasing independence of the planes 
further highlights the magnetic frustration introduced by 
doping. With increasing doping, as the layers become 
decorrelated magnetically, in terms of the superconduc- 
tivity the material becomes more strongly coupled. 

Another explanation for the presence of optic weight at 
low-energies has been outlined in recent theoretical work 
based on a bilayer stripe model by Kruger and Scheidl»i& 
They considered three different stripe structures based 
on the position of the rivers of charge in one layer with 
respect to the other. For each structure they found that 
a finite bilayer coupling would introduce optic spectral 
weight at energies much lower than the optic gap of the 
parent insulator and quite close to the resonance energy 
of the doped system (defined here as the energy at which 
the incommensurate branches of spin excitations meet). 
Such an interpretation could explain the apparent soft- 
ening of the optic mode with increasing doping and the 
presence of a weak resonance in the optic channel.'*'^ We 
note that these calculations are based on an ordered mag- 
netic ground state and therefore do not predict the large 
suppression of spectral weight at low-energies which oc- 
curs for non-zero doping, nor the enhancement of spectral 
weight around the resonance energy. The prediction of 
the bilayer stripe model, that the optic scattering would 
be incommensurate, has been discussed in relation to our 
data in the previous section. 
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FIG. 14: Momentum integrated susceptibility as a function 
of energy transfer for the acoustic mode only. The integra- 
tion for the MAPS data (Ei = 150 meV) is plotted above 35 
meV with filled circles. The open circles are data taken 
previously at Chalk River using a triple-axis spectrometer 
(E/=14.5 meV). The low-energy spectrum is dominated by 
the resonance peak at ~ 33 meV. The high-energy integrated 
susceptibility is nearly constant and is similar to that previ- 
ously measured in the YBCOe.is insulator. 



IV. INTEGRATED INTENSITY AND THE SUM 
RULE 



One of the surprising results from our reactor study 
of the low-energy scattering was that remarkably little 
spectral weight resides at energies less than or equal to 
the resonance energy. We found that the total spectral 
weight made up only about 3% of that expected from the 
total moment sum rule. If the integrated spectral weight 
were distributed equally over all energies, as it is in the 
insulating antiferromagneti^, up to the peak of the dis- 
persion at ^ 250 meV, then we would have expected to 
have found 40 meV/250 meV ~ 16% of the total spec- 
tral weight. The presence of so little spectral weight at 
low-energies necessitates a detailed evaluation of the in- 
tegrated intensity at energies above the resonance energy. 

To get a qualitative idea for the amount of spectral 
weight we observe, and the fraction of the total moment 
which resides around the resonance energy, we consider 
the sum rule for two uncoupled Cu02 layers. Ignoring the 
effect of the chains, and therefore taking only two Cu^+ 
ions per formula unit, the total spectral weight integrated 
over all energy and momentum should be given by the 
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total moment sum rule: 

j(L>Jd\S{q,u;)=2x^S{S+l)g^ (3) 

where J d^q is the momentum space integral over the 
correlated peak. In this equation the factor of 2 comes 
from the fact that we are taking two Cu^+ ions per for- 
mula unit. Here we have assumed that ^(q, a;) has been 
corrected for the bilayer structure factor and the Cu^+ 
anisotropic form factor. We do expect the total moment 
sum rule to be reduced as some of the spins will be in- 
evitably destroyed by the doped holes. Since the doping 
concentration is so small, p ~ 0.1 in the case of Ortho-II 
YBCOg.s, we expect the total moment sum rule to be 
approximately obeyed. Equating the above integral to 
the cross-section for paramagnetic scattering 

1 = ^-^ j duj j d\[n{uj) -f l]x"(q, u:) = (4) 

For S=i and g — 2 this gives a total integral of 1= 
This estimate provides a useful benchmark with which 
we can compare our integrated intensities. 

By integrating the data from to 120 meV in both 
momentum and energy over the optic and acoustic 
modes, we find a total integral of 7r~^ J d^q J du! [n{uj) + 
l](Xac.(q,w) + Xop.(q,^)) ~ 0.3 /x| or about 15 % of the 
total weight predicted by the sum rule. The integral from 
40 to 120 meV then gives ^ 0.23 /Lt^. As a consistency 
check we have integrated the acoustic mode from to 40 
meV and have found a total spectral weight of ^ 0.07 
IJ,^ which agrees reasonably well with the total integral 
of 0.05 measured in our reactor experiment. The dif- 
ference between the reactor and spallation absolute unit 
calibration illustrates that the experimental error in the 
measurements quoted here is about 20% or ^ 0.015 /x^. 

To connect our data with other experiments and the- 
ory, we must correct the absolute measurement for the 
quantum renormalization factor Z^. has been mea- 
sured from the momentum-integrated spectral weight in 
both insulating YBCO and LSCO and agrees with the 
calculated value of ~ 0.5 ASiS Therefore, the actual frac- 
tion of the total moment measured is about 15 % of the 
total moment. In the insulator YBCOg.is, performing 
the same integral up to 120 meV one gets ^ 0.2 fi^/eV.'^ 
It is not surprising that we obtain a small fraction of the 
total moment at these energies as x" increases dramat- 
ically as the dispersion reaches its maximum near the 
antiferromagnetic zone boundary. 

A simple prediction of linear spin-wave theory 
for a square plaquette of Cu^^ spins is that the 
low-temperature momentum-integrated susceptibility, 
x"(a;) = J d^q x"(q, w), is independent of uo in the lin- 
ear spin-wave region as has been confirmed in insulating 



YBCOe.ig^^. As pointed out by Shamoto et at, the 
momentum-integrated weight is inversely proportional to 
the exchange constant44 Therefore, a reduction in the 
exchange constant (or spin-wave velocity) will translate 
into an increase in the momentum integrated spectral 
weight. For the insulator, in the low-temperature or- 
dered phase, x"{^) saturates at a value of ~ 4 fi'^/eV. 
The momentum integral J d^ qx" to) as a function of 
energy transfer for YBCOg.s is plotted in Fig. ^| The 
MAPS data are shown for the range from 35-120 meV and 
our earlier Chalk River reactor data are plotted from 0- 
45 meV. The reactor rather than MAPS data are shown 
below 45 meV because the poorer energy resolution of 
MAPS for Ei=150 meV causes the resonance to appear 
much broader in energy and lower in peak intensity than 
that measured in our reactor data. We emphasize, how- 
ever, that the total integrals for the reactor and time-of- 
flight measurements are the same within error for energy 
transfers 0-45 meV. The two data sets do agree in the 
energy range of 35-45 meV where the magnetic scatter- 
ing is smoother in energy. As illustrated in Fig. ^1 the 
prediction for the spectral weight from linear spin-wave 
theory is violated at low-energies in YBCOe.s where the 
spin fluctuations are dominated by the intense resonance 
peak. 

Because the spectral weight at low-energies has been 
fully discussed in Ref. 0, we will focus on the spectral 
weight above the resonance energy. Fig. E| shows that 
the momentum integral above 40 meV is nearly constant 
as expected for spin excitations with a linear dispersion. 
Its magnitude is consistent with measurements made by 
Fong et aP^ (see Fig. 5). However, our data provide a 
strong indication that x"i'-^) is effectively constant well 
above the resonance energy, unlike previous interpreta- 
tions of a broad second peak located at a high energy of 
60 meV. The high-energy acoustic momentum integral is 
only slightly larger than that (4 ^"^/eV) measured in the 
insulator YBCOg.is by Hayden et al. at 100 meV4i This 
is surprising given that the effective exchange constant, 
and hence velocity, in YBCOg.s has decreased by a large 
factor of 60 %. This would imply (given the same spec- 
tral weight distribution) that the data should approach 
a constant value of 4/0.6 ~ 7-8 ^^/eV. Our data in the 
range 50-100 meV lie in the range 5-7 somewhat 
smaller than expected but within error. A small reduc- 
tion in the renormalization factor may possibly occur 
in the doped system. 

For the total moment, measurements on YBCOg.is 
have given about ~ 0.2 fi^ when integrating up to 120 
meV, whereas we find a slightly larger value of 0.3 /i^. 
Because our measured exchange constant is reduced to 
60 % of that of the insulator we expect our measured to- 
tal spectral weight to be ~ 0.3-0.4 /i^. Our result for the 
total spectral weight is in good agreement with expecta- 
tions given that the accuracy of the intensity calibration 
is about 20 % (as well as the accuracy of the estimated 
effective exchange constant). We emphasize that the to- 
tal integral would be inconsistent, and too small, were 



15 



it not for the fact that the optic mode softens consider- 
ably from the insulator concentration of a;=0.15 to a;=0.5. 
Therefore, the apparent loss in spectral weight displayed 
in the acoustic momentum integral is compensated by 
the presence of optic spectral weight at low energies. For 
Lai.875Bao.i25Cu04, no correction for any change in ex- 
change constant is required since the spin- wave velocity is 
similar to that of the insulator.® The high-energy momen- 
tum integrated spectral weight in Lai.875Bao.i25Cu04 is 
similar to that measured in the insulating monolayer 
compounds. In conclusion, we find that the spectral 
weight at high-energies in YBCOg.s to agree very well 
with that of the insulator when the change in exchange 
constant and the softening of the optic mode is taken into 
account. 

Despite the fact that the resonance carries a large frac- 
tion of the spectral weight at low-energies, it is not sig- 
nificant compared with the total integral extending up to 
higher energies. This lends support to the idea that the 
resonance itself cannot be associated with a pairing boson 
which would require a much larger fraction of the total 
spectral weight to be centered at the resonance energyJS^ 
A low spectral weight around the resonance is predicted 
by the spin-fermion model based on phase space argu- 
ments."*^ 

As an intuitive guide we have carried out linear spin- 
wave calculations based on both an ordered and disor- 
dered static stripe ground state (see Appendix). The 
hole-rich region is modelled as a vacant site across which 
a weak exchange coupling acts to produce antiphase do- 
mains. We expect the model to be most appropriate for 
high frequencies where the excitations see the slow stripe 
dynamics as essentially static. The stripe spin waves may 
provide a partial explanation for the presence of reduced 
spectral weight around the (tt, tt) position. The stripe 
model with small damping (Fig. I15|l shows that signifi- 
cant spectral weight is located throughout the Brillouin 
zone, in contrast to the conventional spin wave model 
whose weight is concentrated around (tt, tt). For large 
damping, which we speculate may occur for YBCOe.s, 
the extra spectral weight would be distributed through- 
out the entire zone, would be unobservable as a well- 
defined correlated peak and so would not contribute to 
our integral for the spectral weight. Although this would 
entail a reduction of weight near (tt, tt), our calculations 
show that the weight located far from the (tt, tt) position 
(and below huj ~ 120 meV) comprises only ~ 10 % of 
the total spectral weight. Therefore, in terms of spectral 
weight located around (tt, tt), the stripe model does not 
predict a substantial change from that of the insulator. 

The dispersion and distribution of spectral weight at 
high energies predicted by the stripe calculation shown 
here (and indeed those presented elsewhere) differ dra- 
matically from the data. Linear spin-wave theory based 
on an ordered stripe ground state predicts bands in the 
high-energy spectrum with substantial gaps of order 0.5 
50 meV below Uui ~ J. Instead of large gaps in the 
energy spectrum we observe the spectral weight to be 



distributed uniformly above the resonance energy. The 
predicted spectrum based on an ordered stripe ground- 
state would also result in an significant anisotropy in the 
high-energy dynamics. This is clearly not the case in 
our measurements as illustrated in Figs. (Sjand^and as 
discussed previously. Since we know the measured high- 
energy spin spectra are continuous in energy as well as 
relatively isotropic (Fig. [SJ we have extended our stripe 
model to include roughened stripe boundaries. This may 
also simulate stripes moving slowly relative to the high- 
frequency spin waves. We have performed calculations 
where the position of the hole-rich region, represented by 
a vacancy, is randomly placed to the left of, to the right 
of, or on the mean position of the boundary in the above 
8-spin-plus-vacancy stripe model. The results given in 
Fig. 1161 show that the disordered stripe model partially 
eliminates the band-gaps, however there still exists a 
large anisotropy not observed in our data. We therefore 
conclude that only the low-energy dynamics are similar 
to that expected from an ordered stripe ground state. 

Although the stripe model correctly predicts incom- 
mensurate modulation below the resonance, we empha- 
size that the model, in particular its spectral weight, is 
not valid at low energies. This is because the stripes are 
dynamic (Ref. rather than static in YBa2Cu306.5, 
because the bilayer coupling is assumed small, and be- 
cause spin destruction by hole doping is ignored. The 
model provides only a qualitative way of interpreting the 
spin map in q and lu at high energy. 

The spectral weight for a ladder is illustrated in the 
recent calculations of Uhrig et al?'^ who calculated the 
triplon excitation spectrum associated with a two-leg 
spin ladder. Two van Hove singularities were found, 
around the resonance energy (or when the incommen- 
surate branches meet) and at the top of the spin en- 
ergy dispersion. It is likely that a good fraction of the 
remaining spectral weight exists at higher energies (the 
high-frequency van Hove singularity) and therefore is not 
included in our total moment integral. The calculations 
of Uhrig et al. also show that the momentum-integrated 
spectral weight is nearly constant between the two van 
Hove singularities. This agrees with our observations. 
However, their calculations predict the absolute spectral 
weight at high-energies to be substantially reduced from 
that in a generic cuprate insulator (like YBCOg.is). This 
is inconsistent with our results and those measured in 
other cuprates. Experiments extending to much higher 
energies are clearly required to resolve the issue of total 
spectral weight. 



V. DISCUSSION 

One of the key observations in this experiment is the 
presence of excitations at high-energy transfers similar 
to the spin- waves observed in the parent insulating com- 
pounds. This is based on the high-energy dispersion 
as well as the absolute momentum integrated spectral 
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FIG. 15: (Color online) The spectral weight calculated for an 
ordered spin structure with domains of eight spins separated 
by a stripe of 5 = states. The dispersions along the a* 
direction is shown in the upper and panel and along the b* 
in the lower panel. The intensity is plotted on a logarithmic 
scale. The value of a and the domain size was chosen to agree 
with our reactor data, as discussed in the text. 
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FIG. 16: (Color online) The spectral weight calculated for 
a disordered stripe state on a logarithmic scale. Instead of 
the hole always occupying the same position, the hole was 
allowed to shift randomly along [100] perpendicular to the 
stripe boundary. The detailed bands found in the ordered 
stripe state are found to be smeared out considerably. 



weight at high-energies. The spin excitations are well- 
defined in wave vector but not in energy. A direct com- 
parison of these resuhs to those of the low-energy excita- 
tions previously studied suggests two distinct regions of 
the excitation spectrum: one at low-energies, less than 33 
meV, and well described by one-dimensional incommen- 
surate fluctuations, and another at high-energy with ex- 
citations characterized by commensurately centered fluc- 
tuations that disperse in a similar manner to spin- waves. 

Based on the measured dispersion in 
Lai.875Bao.i25Cu04" and Lai.84Sro.i6Cu04-^, the 
low-energy dispersion with incommensurate peaks meet- 
ing at higher energies appears to be a common trait of 
the cuprate superconductors. However, the momentum 



integrated acoustic susceptibility as a function of energy 
transfer we measure in YBCOg.s contrasts with that 
measured in Lai.86Sro.i4Cu04i^ In particular no sharp 
well-defined resonance peak in the susceptibility is 
observed in LSCO. As noted in our previous reactor 
study, the spectral weight under the resonance comes 
from the low-energy spectral weight suppressed in the 
superconducting state. Therefore, a combination of the 
a higher energy at which the incommensruate peaks 
meet and a lower superconducting gap may explain 
why the resonance peak is not as well defined in LSCO 
compared with YBCO. 

The occurrence of incommensurate scattering that 
is described by a dynamic stripe model only at low- 
energies suggests that stripes are weakly pinned in the 
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YBCOe+a; system. This result provides an important 
energy scale for theories involving low-energy fluctuating 
stripe phases4£ The evolution with energy from stripe- 
like excitations that are hydrodynamic in nature to spin- 
waves also provides an explanation for the breakdown 
of a;/T scaling observed at high-energies in our previous 
reactor work. We emphasize that even though wc use 
the language of stripes here, any model which results in 
one-dimensional spin density wave fluctuations at low- 
energies would be consistent with our data. 

One way of reconciling these two regions in the excita- 
tion spectrum in terms of fluctuating stripes is that the 
resonance may represent a characteristic energy defined 
by the size of a stripe domain. For spin fluctuations with 
a wavelength greater than the domain size (low-energy 
excitations) the dynamics of the domain walls will be 
important. For high energies where the wavelength is 
shorter than the distance between domain walls the dy- 
namics of the stripes will not be important and the exci- 
tations will look similar to those of the parent insulating 
compound. This physical picture is partially reflected in 
our mean-field calculations (see Appendix) based on a 
disordered stripe ground state. This response is quali- 
tatively similar to that of stripe-liquid theoryii^ The ex- 
periment shows that the spin response then shows stripe- 
like behavior only at low energies. This then implies, in 
the context of stripes, that the rivers of charge are only 
weakly pinned. 

The excitation spectrum of a stripe ground state was 
first calculated by Batista et alw^ and has recently been 
extended by Kruger and Scheidli^ and Carlson et alm^ 
These calculations, as well as our own linear spin- wave 
calculation, show that a stripe ground state cannot ex- 
hibit, in the case of weak stripe coupling and disorder 
(see Appendix), a symmetric cone around the (tt, tt) po- 
sition. Despite the fact that disorder helps to remove 
the gaps, a significant anisotropy still exists in the line 
width at high-energies when scanned along [10] compared 
with [01]. This anisotropy results from the intersection of 
two spin- wave cones emanating from the incommensurate 
positions. This is in contrast to the relatively isotropic 
widths we observe in YBCOg.s. Also, the spin- wave cal- 
culation is not capable of explaining the suppression of 
spectral weight at low-energies, for it ignores the mag- 
netic frustration introduced by doping and as well as- 
sumes an ordered ground state. 

The existence of excitations above the resonance en- 
ergy, similar to those found in the parent insulating 
cuprate compounds, suggests that the resonance may be 
interpreted as a gapped mode which gradually decreases 
in energy with decreasing doping. Such an energy gap 
was predicted in the early theory of Chubukov et alM for 
a quantum disordered ground state. More recent develop- 
ments by Demler et al. suggesting the close promximity 
of a spin density wave (SDW) phase may explain the ex- 
istance of the low-energy incommensurate fluctuations.^^ 
The magnitude of the energy gap is predicted to decrease 
according to a power-law as the Neel state is approached. 



which is qualitatively seen in the YBCOg+a; system. Such 
a response is also predicted in the spin-fermion model^ 
To test this prediction it would be valuable to study the 
high-energy excitations near the boundary of antiferro- 
magnetism and superconductivity. 

Predictions for the high-energy dispersion have been 
made in the random-phase approximation by Brinck- 
mann and Lee'*^ and by Kao et al.^^ Brinckmann and 
Lee predict the scattering to be incommensurate below 
the resonance and to be commensurate at the resonance 
energy, consistent with experiment. The dispersion above 
the resonance is predicted to have a separation between 
+q and —q branches in an (H 0) scan, that increases 
with energy, similar to spin waves 4i In the model by 
Kao et al. the separation in q is constant or slightly de- 
clines with increasing energy (see Fig. 1 of Ref. IT^ . 
unlike our observation of an opening in the q separation 
with energy that reflects a characteristic spin velocity. 
The high-energy line shape predicted in Ref. ^| is nearly 
isotropic and is consistent with the near symmetric ring 
observed in this experiment. 

Despite the several similarities in the high-energy spec- 
trum of the cuprates recently measured and discussed 
in this paper, there are many key differences which are 
currently not resolved. First, despite the fact that we 
observe a symmetric ring of scattering, recent measure- 
ments reported in Ref. (3 show a clear anisotropy around 
the (tt, tt) position. Also, despite the fact we observe 
the momentum integrated spectral weight to be nearly 
frequency independent above the resonance energy, pre- 
vious studies (including Ref. and Ref. IsgI) report a 
significant peak in the integrated spectral weight at high 
energies. Some of these apparent discre pan cies may be 
associated with different dopings as Ref. 0| has suggest- 
ted the existance of a metal to insulator transition in 
YBCOe+a; at 0.07. Nevertheless, further studies are 
required to resolve these discrepancies. 



VI. CONCLUSION 

We have shown that the spin excitations of an under- 
doped YBCO superconductor evolve from stripe modu- 
lations at low-energies to a near symmetric cone of spin- 
waves above the resonance energy. Based on the disper- 
sion we have estimated the in-plane exchange constant 
to be ^ 75 meV, 40% less than the exchange of 125 meV 
in YBCOg.is. In contrast to the insulating compound 
we find that optic weight is present at relatively low- 
energies. We suggest that with increasing doping the bi- 
layer spins are gradually decoupled as the superconduct- 
ing coupling of the carriers increases. We find that the 
total spectral weight is similar to that expected from the 
insulator, when renormalization factors and the reduced 
exchange constant are taken into account. The form of 
the high-energy spectrum cannot be explained with an 
ordered stripe model, nor one with weak disorder. The 
high-energy excitations are similar to spin-waves in the 
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insulator based on both dispersion and spectral weight. 
The data are inconsistent with what is expected from a 
static ladder ground state. We note that the several sim- 
ilarities of the excitation spectrum we measure in Ortho- 
II YBCOg.s with that measured in Lai.875Bao.i25Cu04 
suggests that the spin dynamics in doped Cu02 systems 
may have a universal explanation. 
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VII. APPENDIX: EXCITATIONS OF AN 
ORDERED AND DISORDERED STRIPE 
GROUND STATE 

The qualitative nature of the spin excitations above 
and below the resonance energy are very different. The 
low-energy excitations are incommensurate with peaks 
displaced only along the [100] direction, in qualitative 
agreement with the spectrum predicted for a stripe 
ground state. The high energy excitations are isotropic 
in the direction of q, within experimental error, and are 
qualitatively similar to spin-waves originating from the 
(tt, tt) position. Given the different behavior at high 
and low-energy energies it is important to determine if 
the two regimes can be understood in terms of the same 
ground state. In this section we describe linear spin- wave 
calculations for an ordered stripe ground state as well as 
numerical simulations describing the effects of disorder 
on the spectrum. We find the results of the calculation 
are not consistent with experiment. 

In our calculations, which parallel those of Batista et 
alm^, Kruger and ScheidU^ti^, and Carlson et alM: we 
have numerically calculated the excitation spectrum of 
an ordered stripe ground state with one-dimensional do- 
mains of 8 spins separated by a vacant spin-site. The 
vacant site is taken to model the region of large hole 
density separating antiphase domains. The size of this 
domain structure is expected from the known position'^ 
of the low-energy incommensurate modulation in Ortho- 
II YBCOg.s. We have used antiferromagnetic Heisenberg 
exchange between the spins and have varied the strength 
of the weak coupling between domains across the vacancy 
at the charge stripe, (K), relative to the strong couphng 
within a spin domain, (J), so as to define the key ratio 
K/ J = a. Having fixed the domain size from the position 
of the low-energy incommensurate peaks the parameter 



a determines the position at which the low-energy in- 
commensurate rods meet, or where the resonance energy 
occurs in experiment. The weak antiferromagnetic ex- 
change, e.g. a = 0.075, between domains couples spins 
two sites apart and so leads to antiphase domains. We 
have determined the relative cross-sections as well as the 
energies so as to compare with experiment. 

The energy spectrum of an ordered spin ground state 
can be determine from dynamical mean-field theory by 
solving the set of coupled, first-order equations given by 

dt{S{l,t)) ^h{l,t)x {S{l,t)). (5) 

Here, h{l,t) is the molecular field acting on the spin 
S{l,t) caused by its nearest neighbors, 

h(/,i) = ^Ja,Z + nO(Sa + (6) 

By assuming a particular ground state structure (S°(^)) 
and expanding the above equations linearly in deviations 
about this ground state {5{S{l,t))) the following set of 
linear equations is obtained 

-zdtS+{l,t) = (7) 

3 

or, 

-idtS+{t) = AS+{t). (8) 

In this equation S+{l,t) = S{S^{l,t)) + iS{Sy{l,t)), 
is the molecular field with (S(Z + ni,t)) set to the ground 
state values ((S*'(Z))), and S{l,j) is equal to one when 
I = j and zero otherwise. By diagonalizing the matrix A 
the energy spectrum can be found for any ground state 
and spin- coupling. 

To obtain the spectral weight as a function 
of momentum and energy transfer we have solved 
for the Green's function G{S+{l,t), S' {j,0)) = 
-ie{t){[S+{l,t),S-{j,0)]}, and then taken the imagi- 
nary part to get the neutron scattering cross-section. 
Solving the equations of motion for the Green's function, 

J2m, k)u; + A{1, k))G{S+{k, u;),S-{j, 0)) = (9) 

k 

{[S+{l,0),S-{j,0)]), 

the neutron cross-section was obtained by setting 
Lu to Lu + ie, inverting, and summing all possible 
G{S'^{l,uj),S^{j,0)). The parameter e represents some 
form of damping. 

The results for the eigenvalues and the spectral weight 
are shown in Fig. 1151 for a domain size of 8 spins and a 
ratio of inter- to intra-domain coupling of a=0.075. Tak- 
ing J=125 meV, we have chosen the domain size to give 
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the correct position for the low-energy incommensurate 
peaks as measured in our reactor data. The value of a 
was chosen such that the top of the lowest branch of the 
dispersion coincided with the position of the measured 
resonance energy. ■^■-'^^ We note that such a detailed micro- 
scopic model for the resonance has not been conclusively 
proven to date. 

The spectrum displays eight branches, as expected 
given the ground state chosen in this calculation. The 
energy axis is given in units of the intra-stripe coupling 
(J). It is worth noting that the spectral weight is con- 
centrated roughly in the region where excitations would 
be present in the insulating compound. Thus, we see 
in Fig. El that the weight lies close to where a single 
conventional spin-wave band would occur as one moves 
from one band to another of the multi-branch dispersion 
of the stripe model. The concentration of spectral weight 
around individual branches separated by a finite gap is 
not observed in the experiment. It is satisfying, how- 
ever, that the spin waves of the stripe model do show 
one-dimensional anisotropy similar to observation: along 
[0.5, K] a gap occurs at the resonance energy with little 
weight at lower energy, while along [H, 0], perpendicular 
to the stripe, the energy falls to zero and the intensity 
grows on approach to the incommensurate momenta. 

Based on our reactor studies of the low-energy spec- 
trum of Ortho-II YBCOg.s, we expect the spin response 
to be highly damped. In our previous analysis of the 
low-energy scattering we estimated a dynamic correla- 
tion length K 20 A and a spin-wave velocity of he ~ 
300 meV • A giving a large energy broadening of 2T ~ 30 
meV. In terms of the stripe picture this could be taken 
as evidence for the presence of a large amount of disorder 
along the stripe of 5=0 states. In an attempt to simu- 
late this we have arranged the holes on a finite lattice to 
have a random and normal distribution (with standard 
deviation a=l cell or one Cu^"*" site) around the ordered 
position such that the hole can be shifted up two unit 
cells along the [1, 0] direction. Using periodic boundary 
conditions we have then calculated the spin spectrum 
based upon a 1024 spin lattice ground state(32 by 32 



spins) using the same formulae presented above. The re- 
sults were not found to change significantly for smaller 
lattices. The results presented in Fig. ^jare averaged 
over seven different iterations of the 32 x 32 random- hole 
structure. 

The results for such a configuration is presented in Fig. 
Eland shows that the spectral weight can be smeared out 
considerably by such a disordered ground state. This cal- 
culation also shows that many of the features displayed 
by the ordered system displaced far from the (tt, tt) po- 
sition are smeared out. Also, the spectral weight again 
appears as a correlated peak around the (Q, w) positions 
expected in the absence of stripes and the gapped struc- 
ture calculated for the ordered system is now smeared 
out in energy and momentum. Because the model as- 
sumes an ordered ground state it is not able to repro- 
duce the suppression of the low-energy spin fluctuations, 
which we have previously found to be well described by 
overdamped spin-waves. Our formulation of a disordered 
stripe structure is expected to be more appropriate for 
higher-energy excitations (above the resonance) as the 
energy transfer is greater than the broadening of 2 F 
~ 30 meV expected from our reactor data. We note 
that the high-energy response found in our disordered 
stripe model is consistent with the increase with energy 
of the observed width in q reported for commensurate 
Lai.875Bao.i25Cu04, although our data and calculations 
show less anisotropy at high energy than the calculations 
for the spin ladder modeL- 

Even though the low-energy scattering below the reso- 
nance may be interpreted in terms of stripes, our numer- 
ical calculations show that it is difficult to interpret the 
high-energy scattering in terms of a stripe ground state. 
Our results for an ordered ground state (with no disor- 
der) show large band gaps clearly not observed in the 
data. The introduction of disorder does the band gaps 
but does not remove a large anisotropy, particularly in 
the linewidth. Our stripe model fails to account for the 
peak in the spectral weight at around 55 meV. The model 
does not account for the asymmetric resonance spectral 
profile found in our previous reactor worki^ 
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